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1.  INTRODUCTION 

1.1  General  Introduction 

The  research  program  for  Contract  Number  F 1 9628-82-K-0032  had 
both  scientific  and  technical  objectives.  The  scientific 
objectives  included  the  measurement  of  the  polarizations  and  wave 
characteristics  of  Pi2s.  These  data  were  then  to  be  used  to 
develop  a  model  of  the  magnetic  fields  and  currents  associated 
with  Pi2s.  This  model  was  to  be  used  in  evaluating  methods  for 
detecting  Pi2s  with  ground  based  magnetometer  arrays.  The 
technical  objectives  included  the  design  of  algorithms  for  the 
estimation  of  the  wave  and  polarization  characteristics  of  Pi2s, 
the  programming  of  these  algorithms,  and  the  writing  of 
documentation  to  accompany  the  computer  programs 

The  work  under  this  contract  falls  in  four  specific  areas 
including : 

A.  Pi2  Experimental  Studies. 

These  studies  used  data  from  ground  based  magnetometer 
arrays,  including  the  AFGL  (Air  Force  Geophysics 
Laboratory  ),  the  Alberta  (University  of  Alberta),  and 
the  IMS  (International  Magnetosphere  Study)  arrays.  The 
studies  were  designed  to  measure  the  polarizations, 
phase  velocities,  group  velocities,  frequencies,  and 
wave  envelopes  of  the  Pi2s. 
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B.  Pi2  Models 

These  models  are  based  on  field  aligned  and  ionospheric 
currents  and  were  used  to  model  the  magnetic  fields  of 
Pi2s  at  ground  level. 

C .  Multivariate  Time  Series  Studies 

These  studies  used  the  methods  of  multivariate  analysis 
to  develop  analytic  forms  for  the  estimation  of  the 
parameters  of  waves,  and  in  particular  Pi2s. 
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1.2  An  introduction  to  Pi2s 

Pi2  pulsations  are  oscillations  of  the  geomagnetic  field  which 
occur  in  conjunction  with  intensifications  of  the  magnetospher ic 
currents  associated  with  polar  magnetic  substorms.  The  Pi2s  are 
seen  as  quasi  sinusoidal  oscillations  on  ground  based 
magnetometers.  The  oscillations  typically  have  frequencies  of 
about  10  mHz  (0.01  Hz),  and  amplitudes  of  up  to  several  nT  at  mid 
latitudes  and  tens  of  nT  at  high  latitudes.  A  pulsation  train  can 
last  for  one  cycle  or  more. 

These  pulsations  are  probably  the  magnetic  signatures  of 
transient  changes  in  the  ionospheric  and  field  aligned  currents 
(fac)  associated  with  the  onset  of  the  substorm  expansive  phase. 
Observations  which  support  this  possibility,  are  those  showing 
the  association  of  Pi  2's  with  the  initial  brightening  of  auroral 
arcs  and  precipitation  of  electrons  at  the  beginning  of  the 
substorm  expansive  phase  [Troitskaya  and  Gul'elmi,  1967;  Pytte 
and  Trefall,  1972;  Samson  1982],  and  the  close  correlation  of  the 
polarization  -  pattern  of  Pi  2's  with  the  magnetic  fields  from 
the  substorm  current  "wedge"  [Rostoker,  1967;  Bjornsson  et  al., 
1971;  Baransky  et  al.,  1980;  Lester  et  al.,  1983,  1984;  Samson 
and  Harrold,  1983].  (See  McPherron  et  al.  [1973]  for  a 
description  of  the  substorm  current  "wedge".) 

Samson  and  Harrold  [1983]  showed  that  all  the  various  ground 
based  observations  of  the  polarizations  could  be  fitted  to  a 
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consistent  pattern,  with  counterclockwise  (cc)  polarizations  at 
mid-latitudes,  and  a  more  complicated  pattern  near  the  region  of 
the  onset  of  substorm  fac  (see  Fig.  7  in  Samson  and  Harrold) . 
Studies  conducted  under  this  contract  [Samson,  1985]  have  shown 
that  the  observed  pattern  of  polarizations  is  compatible  with 
westward  moving  fac  and  ionospheric  current  systems  which  are 
associated  with  a  brightening  arc. 

The  midlatitude  wave-numbers  and  phase  velocities  of  Pi2's 
[Mier-Jedrze jowicz  and  Southwood,  1979,  1981;  Baransky  et  al., 
1980;  Lester  et  al.,  1983,  1984]  indicate  that  most  Pi2’s  show 
westward  propagation,  with  phase  velocities  near  2  degrees 
longitude  /sec  (  longitudinal  wavelengths  near  200  degrees,  and 
mean  frequencies  near  0.01  Hz)  at  least  near  the  longitudinal 
center  of  the  substorm  current  wedge. 
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2.  WORK  COMPLETED 


2.1  A.  Pi2  experimental  studies 

The  experimental  work  used  data  from  the  AFGL,  Alberta,  and 
IMS  arrays  (Fig.1)  to  study  the  polarizations,  phase  velocities, 
group  velocities,  and  characteristic  time  constants  of  Pi2s. 

The  Pi2s  appear  to  be  correlated  with  the  brightening  of 
auroral  arcs  [Samson,  1905].  Often  this  arc  is  equatorward  of  the 
Harang  Discontinuity.  The  Pi2s  show  their  maximum  intensities 
near  the  brightened  arc  and  the  substorm  enhanced  auroral 
electrojets  (Fig. 2).  Although  the  polarization  patterns  of 
individual  events  can  be  quite  complicated,  these  patterns  do 
appear  to  be  consistent  from  event  to  event  (Fig.  3,  see  also 
Fig.  7  in  Samson  and  Harrold  [1983]).  Near  the  longitudinal 
center  of  the  substorm  current  wedge,  the  polarizations  range 
from  counterclockwise  (cc)  (viewed  down  on  the  H-D  plane)  at 
midlatitudes,  to  clockwise  (cw)  just  equatorward  of  the  intensity 
maximum,  to  cc  just  poleward  of  the  intensity  maximum.  Our  model 
studies  have  shown  that  these  polarizations  are  compatible  with  a 
westward  moving  region  of  transient  fac  and  ionospheric  current 
associated  with  a  brightening  arc  [Samson,  1985]. 

Analyses  of  data  from  the  midlatitude  AFGL  stations  have 
yielded  excellent  measurements  of  the  characteristic  rise  and 
fall  times,  frequencies,  phase  velocities,  and  group  velocities 
of  Pi2s  [Samson  et  al.,  1985].  These  measurements  were  based  on  a 
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polarized  wavelet  model  of  the  vector  time  series  of  the  Pi2  (see 
the  discussion  below  or  Samson  et  al.[l985].).  Measured  group 
delays  over  the  AFGL  array  were  often  less  than  5-10  sec, 
indicating  that  Pi2s  should  give  good  timings  of  substorm  onsets. 
The  rise  and  fall  times  for  the  envelopes  of  the  Pi2  pulsation 
trains  were  typically  60-100  sec. (Fig.  4).  These  short  time 
intervals  suggest  that  the  Pi2s  are  not  associated  with 
hydromagnetic  resonances.  However,  the  rise  and  fall  times  are 
quite  compatible  with  transient  reflections  of  fac  over  a 
brightening  arc  [Samson,  1985]. 

Near  the  center  of  the  substorm  current  wedge  (see  Samson  and 
Harrold  [1983]  for  a  description  of  substorm  centered 
coordinates),  most  midlatitude  Pi2s  have  westward  phase 
velocities  with  typical  values  of  1  to  4  degrees  longitude/  sec. 
(Fig. 5).  There  are  some  indications  of  eastward  propagation  more 
than  30°  east  of  the  center  (see  also  Lester  et  al.  [1984]). 
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2.2  B .  Pi2  models 


Due  in  part  to  the  inherent  difficulties  in  a  self  consistent 
theoretical  description  of  transient  hydromagnetic  processes  in 
the  magnetosphere,  we  chose  an  heuristic  model  for  Pi2s.  This 
model  was  based  on  fac  and  ionospheric  currents  compatible  with 
reflection  of  an  Alfven  wave  from  a  conducting  strip  in  the 
auroral  ionosphere.  This  conducting  strip  is  oriented  east-west, 
and  is  associated  with  a  brightening  arc.  The  model,  shown  in 
Fig.  6,  has  fac  at  the  edge  of  the  conducting  strip,  and 
ionospheric  currents  closing  the  circuit.  The  model  is  compatible 
with  an  incident  Alfven  wave  with  an  eguatorward  directed 
electric  field.  To  maintain  a  curl  free  electric  fie’d,  this 
configuration  requires  the  longitudinal  scale  size  of  the  Pi2  to 
be  much  larger  than  the  latitudinal  scale  size.  Observations 
indicate  that  this  should  be  the  case  since  longitudinal 
wavelengths  of  Pi2s  are  often  near  180°  (see  Samson  et  al. 

[1985]  and  references  therein). 

More  details  on  the  formation  of  fac  due  to  reflections  from  a 
nonuniform  ionosphere  can  be  found  in  Samson  [1985].  Southwood 
and  Hughes  [1985]  have  also  considered  a  Pi2  model  which  takes 
into  account  reflections  fron  a  nonuniform  ionosphere. 
Unfortunately  their  model  assumes  an  infinite  latitudinal  extent, 
and  longitudinal  wavelengths  which  are  comparible  to  the 
latitudinal  width  of  the  enhanced  conductivity  band  associated 
with  an  arc.  These  assumptions  are  not  compatible  with  the 
observations  mentioned  above. 
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To  run  the  model  simulations,  we  have  used  a  computer  program 
which  computes  the  magnetic  fields  associated  with  the  fac  and 
ionospheric  currents  in  the  model  depicted  in  Fig.  6.  A  dipolar 
geometry  is  used  for  the  fac.  We  assumed  westward  propagation, 
and  consequently,  the  current  system  must  move  westward  at  a 
velocity  appropriate  for  the  measured  phase  velocities.  The  test 
for  the  validity  of  the  model  was  based  on  the  prediction  of  the 
observed  polarization  pattern  of  the  Pi2s.  The  predicted  pattern 
is  shown  in  Fig.  7  and  compares  favorably  with  the  pattern  in 
Fig.  3  and  in  Fig.  7  of  Samson  and  Harrold  [1983]. 
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2.3  C.  Multivariate  time  series 

Methods  of  multivariate  time  series  analysis  were  required  in 
order  to  determine  the  characteristic  polarization  patterns  of 
Pi2s,  to  measure  the  group  and  phase  velocities,  and  to  develop 
data-adaptive  filters  to  enhance  Pi2  waveforms  in  noisy  data.  We 
have  also  made  some  preliminary  progress  in  using  multivariate 
detection  theory  for  the  detection  of  Pi2  wavetrains.  Almost  all 
our  theoretical  work  has  been  based  on  the  concept  of  pure  states 
in  multivariate  processes  [Samson,  1983a].  Numerous  tests  have 
shown  that  Pi2s  fit  this  model  very  well  [Samson,  1985;  Samson  et 
al.,  1985].  We  have  also  developed  theoretical  expressions  for 
reducing  sample  biases  in  the  estimators  of  the  pure  states, 
based  on  the  work  of  Samson  [1983b].  These  revised  estimators 
have  been  incorporated  into  the  FORTRAN  programs  written  under 
this  contract. 

The  specific  model  which  we  have  found  suitable  for  Pi2s  has 
the  mathematical  form 

w(x,t)  =  1/2[a(x,t)£(x)exp(-i27rf  ot)  +  c.c.]  +  e{x,t)  [1] 

where  w{x,t)  is  the  vector  wavelet  of  the  Pi2  at  position  x  and 
time  t,  a(x,t)  is  the  amplitude,  g(x)  is  the  complex  polarization 
vector  [Samson,  T983a],  fo  is  the  frequency,  and  e(x,t)  is  the 
noise  vector  [Samson  et  al.,  1985],  The  estimators  of  the 
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amplitude,  a(t),  phase,  0(t),  and  wave-frequency 
respectively 


aj(t)  =  (r^j{t)  +  q^(t))’'^  [2] 


t  2 

0j(t)  =  (t2-t,  +  1)'’  I  tan*  Uq  j  (t)/r  j  (t)  ] , 

t=t, 


t  2 

f 0  =  {t2-t  1  +  1  )  *  ’  Z  f  (t)  , 

t  =  t, 


and 

2fff(t)  =  a j *  * ( t ) [ r j ( t )q j (t+1 )-q j { t ) r j ( t+ 1 ) ] 

(an 

where 

N/2 

r(t)  =  N*’Z  c(k)b(k)2exp(27rikt/N)  +  c.c., 
k=0 


N/2 

g(t)  =  N*’Z  c{k)b(k)gexp(2irikt/N)  +  c.c., 
k=0 


2  =  E*  ’ ^u, 

b(k)  =  2* [E* ’ -S* ’ ]z ( k )/(2’E* ’p) , 


N-1 

z(k)  =  Z  d(t)exp(-2jrikt/N) , 
t  =  0 
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and 

N/2 

S  =  Z  c(k)z(k)2‘(k) . 
k=0 

[11] 

In  the  above  equations,  E  {£=£„)  is  the  noise  spectral  matrix,  • 
u  is  the  eigenvector  corresponding  to  the  largest  eigenvalue  of 
E' ’ *S  E' ’ ^ ,  and  c(k)  is  the  spectral  window.  The  above 
estimators  were  used  by  Samson  et  al.  [1985]  in  their  analysis  of 
Pi2s.  An  illustration  of  the  use  of  the  above  estimators  in  the 
analysis  of  mid  latitude  Pi2s  is  given  in  Fig.  8,9.  Figure  8 
shows  Pi2  pulsations  recorded  by  the  AFGL  magnetometers,  after 
application  of  a  3  mHz  highpass  filter.  Figure  9  shows  the 
estimated  envelopes  a(t),  and  the  frequency  f(t). 
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2.4  D.  Computer  programs  and  documentation 


Under  the  terms  of  this  contract,  four  programs  have  been 
written  in  FORTRAN,  documented,  and  sent  to  AFGL.  These  programs 
are  POLFILTl,  P0LFILT2,  QUADFILT,  and  LHFILT.  All  of  these 
programs  were  designed  to  aid  in  the  estimation  of  the  parameters 
of  propagating  waves,  and  in  particular  Pi2s. 

POLFILTl  consists  of  a  main  program  and  subroutines  for  the 
adaptive  filtering  of  multichannel  time  series.  The  filters 
extract  the  polarized  or  pure  state  by  modulating  the  spectral 
representation  with  a  function  based  on  the  degree  of 
polarization  of  the  multichannel  data  [Samson,  1983a].  These 
filters  also  have  an  adjustment  for  lagged  arrivals  in  the 
wavelets,  whereas  the  more  sophisticated  P0LFILT2  does  not. 
POLFILTl  has  a  number  of  biases  in  its  estimators  and  is  designed 
largely  to  act  as  a  rapid,  inexpensive  filter  for  the  initial 
processing  of  the  data.  For  more  details  read  the  manual 
DOCUMENTATION  FOR  POLFILTl. 

POLFILT2  is  a  more  ambitious  program,  and  uses  a  number  of 
optimal  estimators  to  determine  the  parameters  of  pure  states  or 
polarized  waves.  This  program  can  be  used  to  filter  data  to 
extract  the  time  series  of  the  pure  states,  and  to  estimate  the 
parameters  of  polarized  waves  (eg.  the  parameters  of  the 
polarization  ellipse).  Because  it  is  a  much  more  sophisticated 
program,  P0LFILT2  tends  to  be  much  more  expensive  to  run,  and 
requires  a  much  larger  number  of  input  parameters  than  does 
POLFILTl.  Further  information  on  P0LFILT2  can  be  found  in  the 


manual  DOCUMENTATION  FOR  P0LFILT2. 
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QUADFILT  is  a  FORTRAN  program  for  estimating  the  amplitudes, 
frequencies  and  phases  of  polarized  wave  packets  (see  equations 
[1]  to  [10]).  This  program  uses  a  quadrature  filter  in  the 
frequency  domain  to  determine  the  Hilbert  transform  of  the 
multichannel  data  [Goodman,  1960]  (see  Equation  [7]).  This 
program  has  been  used  extensively  in  the  estimation  of  the  phase 
and  group  velocities  of  Pi2s  [Samson  et  al.,  1985].  For  more 
details  on  this  program  please  read  DOCUMENTATION  FOR  QUADFILT. 

The  fourth  program  written  for  AFGL  is  LHFILT.  This  program  is 
a  filter  to  low  pass,  high  pass,  or  band  pass  selected  time 
sequences.  We  have  used  this  filter  in  detrending  our  data  before 
using  POLFILTI,  POLFILT2,  or  QUADFILT.  Details  on  this  filter  are 
given  in  DOCUMENTATION  FOR  LHFILT. 
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occuring  on  Day  298,  1979  (after  Samson  [1985]).  The  intensity 


-ure  3  Polarization  in  the  H-D  plane  for  the  event  occuring  on  Day  298,  1979.  Clockwise  polarization 
depicted  by  a  dark  ellipse  with  an  arrow.  The  others  are  counterclockwise. 
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Figure  4.  Rise  times  for  mid-latitude  Pi2s  (after  Samson  et  al.  [1985]).  The  rise 
time  is  from  the  l/e*maximum  to  the  maximum. 


The  predicted  polarizations  in  the  H-D  plane  using  the  model  in  Figure  6. 
indicate  the  ratio  of  the  minor  axis  to  the  major  axis  of  the  polarization 
Negative  values  indicate  cw  polarization. 


Filter  1 
H  Component 


Filter  1 
D  Component 


UNIVERSAL  TIME 


UNIVERSAL  TIME 


Figure  8.  A  P12  recorded  at  the  AFGL  array.  The  time  series  have  been 
detrended  with  LHFILT,  and  a  highpass  cutoff  of  3mHz. 
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Figure  9.  The  estimated  envelope,  a(t)  (Equation  [2]),  and  frequency,  f(t) 
(Eqtiation  [5])  for  the  Pi2  in  Figure  8. 


